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Abstract — In this paper we address the actuator/sensor 
allocation problem for linear time invariant (LTI) systems. 
Given the structure of an autonomous linear dynamical system, 
the goal is to design the structure of the input matrix (commonly 
denoted by B) such that the system is structurally controllable 
with the restriction that each input be dedicated, i.e., it can 
only control directly a single state variable. We provide a 
methodology that addresses this design question: specifically, we 
determine the minimum number of dedicated inputs required 
to ensure such structural controllability, and characterize, and 
characterizes all (when not unique) possible configurations of 
the minimal input matrix B. Furthermore, we show that the 
proposed solution methodology incurs polynomial complexity 
in the number of state variables. By duality, the solution 
methodology may be readily extended to the structural design of 
the corresponding minimal output matrix (commonly denoted 
by C) that ensures structural observability. 

I. INTRODUCTION 

This paper is motivated by the dearth of scalable tech- 
niques for the analysis and synthesis of different large-scale 
complex systems, notably ones which tackle design and deci- 
sion making in a single framework. Examples include power 
systems, public or business organizations, large manufactur- 
ing systems, wireless control systems, biological complex 
networks, and formation control, to name a few. Focusing 
on the last case, consider, for instance, the synchronization 
problem in vehicular formations: Given a communication 
topology for inter-vehicle information exchange and the 
individual vehicle (agent) models, we are often interested 
in addressing the following questions: 

• What is the smallest subset of agents (and specifically 
which ones), that need a dedicated input (i.e., a control 
that directly affects a single state variable), such that 
the system is controllable? 

• Similarly, what is the smallest subset of agents (and 
specifically which ones) that need to be equipped with 
a dedicated output (i.e., an output that measures directly 
a single state variable), such that the entire network state 
may be estimated? 

Referring to the vehicular-formation scenario, the different 
agents that require dedicated controls to achieve system 
controllability, are those that play the role of leaders. Under 
infrastructure and operational constraints, identifying the 
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smallest subset of such agents clearly maximizes the effi- 
ciency of the system. The concerns posed above go beyond 
the vehicular-formation example and are applicable to wider 
classes of large-scale multi-agent scenarios. 

To address these problems, we will resort to structural 
systems theory [1], in which the main idea is to reformulate 
and study of an equivalent class of systems for which 
system-theoretic properties are investigated based on only 
the sparsity pattern (i.e., the location of zeroes and non- 
zeroes) of the state space representation matrices. Such an 
approach is particularly helpful when dealing with systems 
parameter uncertainties. Analysis using structural systems 
provides system-theoretic guarantees that hold for almost 
all values of the parameters, except for a manifold of zero 
Lebesgue measure [2], which may further be characterized 
algebraically [3]. Properties such as controllability and ob- 
servability are referred to as structural controllability and 
structural observability in this framework, as they hold in 
general, i.e, for almost all non-zero entries in the state 
space representation. With this, our design objective may be 
precisely formulated as follows: 

Problem Statement 

Given 

x = Ax, (1) 

where A represents the structural pattern of A (i.e., the 
locations of zeroes and non-zeroes only), our goal is 
V\ Design B (i.e., find the structural pattern) with a min- 
imum number of dedicated inputs, such that (A, B) is 
structurally controllable. Stated formally, characterize 
(all) B € R nx P such that 

^j'e{i,--,j>} ,n} Bij 0, 

in other words 

p 

x — Ax + b i:j Uj, ij G {1, • • • , n}, (2) 
j'=l 

where 6j. represents the ij-th canonical vector and 
iij € K represents the j-th control, such that system 
|2]) is structurally controllable, p G N, and there exists 
no other p' < p that satisfies the previous requirement. 
Solution of V\ also addresses the corresponding optimal 
(minimal placement) structural observability output matrix 

'A pair (Ao,Bq) is structurally controllable if and only if Ve > 0, 
there exist completely controllable pairs (Ai,B±) and (A2,B2) with the 
same structure (location of zeroes and non-zeroes) as (Ao,Bq), such that 
\\Ai - A \\ > e, - B |j > e, and \\A 2 - A \\ < e, \\B 2 - B \\ < e. 



design problem by invoking the duality between estimation 
and control in LTI systems. o 

The literature on structured systems theory is extensive; 
see [2,4-6] for earlier work, see also [1] for a recent survey. 
For applications to optimal sensor and actuator placements, 
the reader may refer to [7] and references therein; however, 
these approaches mostly lead to combinatorial implemen- 
tation complexity in the number of state vertices (agents), 
or are often based on simplified heuristic -based reductions 
of the optimal design problems. Systematic approaches to 
structured systems based design were investigated recently 
in the context of different application scenarios, see, for 
example, [8-13]; for instance, in network estimation, as in 
[8,11] , where strategies for output (sensor) placement are 
provided, ensuring only sufficient (but not necessarily mini- 
mal) conditions for structural observability, whereas in [8,13] 
applications to power system state estimation are explored. 
From the structural observability viewpoint, as a key contrast 
to the above approaches, we study the constrained output 
placement problem, specifically, in which the sensors are 
dedicated, in that, they may only measure a single state 
variable. The formulation that is closest to our setup in 
terms of minimal actuator placement arises from biological 
complex networks [13], where the concept of a driving vertex 
is introduced, i.e., a state vertex through which a subset of 
the graph vertices can be controlled. One of the problems 
addressed in that work may be stated as follows: What 
is the minimum number of driving vertices that make the 
entire network controllable? However, the results in [13] 
hold only for the case in which the system graph is strongly 
connected. In contrast, in this paper, in addition to providing 
the minimum number of dedicated inputs (driving vertices) 
for generic system matrices (digraphs), we also characterize 
the set of all possible minimal feasible configurations, i.e., 
allocations with the minimum number of dedicated inputs 
that ensure structural controllability. 

To summarize, in terms of the dedicated input placement 
problem, the main contributions of this paper include: 1) the 
identification of the minimum number of dedicated inputs 
(to ensure structural controllability); 2) characterization of 
all such feasible minimal dedicated input placement config- 
urations (i.e., the design of B up to column permutations); 
3) algorithmic generation of a minimal feasible input con- 
figuration in polynomial complexity (in the number of state 
variables). 

The rest of the this paper is organized as follows. Section|H| 
reviews some concepts and introduces results (some of them 
new) in structural systems theory and establish their relations 
to graph-theoretic constructs. Subsequently, in Section [TIT] 
we present the main technical results (proof outlines being 
relegated to the appendices), followed by an illustrative 
example in Section [TV] Conclusions are presented in Section 
[V| Finally, Section |V| concludes the paper and discusses 
avenues for further research. 

II. PRELIMINARIES AND TERMINOLOGY 

In this section we recall some classical concepts in struc- 
tural systems, introduced in [4]. 



Given a dynamical system ([TJ, an efficient approach to the 
analysis of its structural properties is to associate it with a 
directed graph (digraph) T> = (V, E), in which V denotes a 
set of vertices and E represents a set of edges, such that, an 
edge (vj, Vi) is directed from vertex vj to vertex Vi. Denote 
by X = {xi, ■ ■ ■ , x n } and U = {ui, ■ ■ ■ , u p } the set of state 
vertices and input vertices, respectively. Denote by £ x ,x = 
{(xi,Xj) : [A]ij ^ 0} and E u .x = {{uj,Xi) : [B]^ ^ 0}, 
to define V(A) = (X, 8 X , X ) and V(A, B) = (XUU, £ x ,xU 
E u ,x). A digraph V s = (V S ,E S ) with V S CV and E s C E 
is called a subgraph of T>. If = V, T> s is said to span 
V. A sequence of edges {(vi,v 2 ), (v 2 ,v 3 ), ■■■ , (v k -i,v k )}, 
in which all the vertices are distinct, is called an elementary 
path from V\ to v k . When v k coincides with Vi, the sequence 
is called a cycle. 

In addition, we will require the following graph theoretic 
notions [14]: A digraph T> is said to be strongly connected 
if there exists a directed path between any two pairs of ver- 
tices. A strongly connected component (SCC) is a maximal 
subgraph T>s = (Vs, Es) of T> such that for every v, w € Vs 
there exists a path from wtotu and from w to v. Note that, 
an SCC may have several paths between two vertices and 
the path from v to w may comprise some vertices not in 
the path from w to v. Visualizing each SCC as a virtual 
node (or supernode), one may generate a directed acyclic 
graph (DAG), in which each node corresponds to a single 
SCC and a directed edge exists between two SCCs iff there 
exists a directed edge connecting the corresponding SCCs in 
the original digraph. The DAG associated with T> = (V, E) 
may be efficiently generated in 0(|V| + |.E|) [14], where |V| 
and \E\ denote the number of vertices in V and the number 
of edges in E, respectively. The SCCs in a DAG can be 
characterized as follows 

Definition 1: An SCC is said to be linked if it has at least 
one incoming/outgoing edge from another SCC. In particular, 
an SCC is non-top linked if it has no incoming edges to its 
vertices from the vertices of another SCC and non bottom 
linked if it has no outgoing edges to another SCC. □ 

For any two vertex sets Si , S'2 C V, we define the bipartite 
graph B(Si, S 2 , E Slt s 2 ) associated with D = (V,E), to be 
a directed graph (bipartite), whose vertex set is given by 
Si U S 2 and the edge set E Su s 2 bv E s u s 2 = {(si,s 2 ) € 
E : si e Si,s 2 e S 2 }. 

Given B(Si, S 2 , Es lt s a )t a matching M corresponds to 
a subset of edges in Es ll s 2 that do not share vertices, i.e., 
given edges e = (si,s 2 ) and e' = (s[,s' 2 ) with si,s[ e Si 
and s 2 , s' 2 € S 2 , e, e! € M only if s\ ^ and s 2 ^ s' 2 . A 
maximum matching M* may then be defined as a matching 
M that has the largest number of edges among all possible 
matchings. The maximum matching problem may be solved 
efficiently in 0{^/\S x U S 2 \\E SuSi |) [14]. Vertices in S x 
and S 2 are matched vertices if they belong to an edge in 
the maximum matching M* , otherwise, we designate the 
vertices as unmatched vertices. If there are no unmatched 
vertices, we say that we have a perfect match. It is to be 
noted that a maximum matching M* may not be unique. 

For ease of referencing, in the sequel, the term right- 



unmatched vertices (w.r.t. B(Si, S2, -Esi,s 2 ) an d a maximum 
matching Af*) will refer to only those vertices in S2 that do 
not belong to a matched edge in Af*. 

A. Structural Systems 

Given a digraph, we further define the following special 
subgraphs: 

- Consider T>(A) : State Stem - An elementary path composed 
exclusively by state vertices, or a single state vertex. 

- Consider T>(A, B) : Input Stem - An elementary path 
composed of an input vertex (the root) linked to the root 
of a state stem. 

- Consider T>(A) : State Cactus - Defined recursively as 
follows: A state stem is a state cactus. A state cactus 
connected to a cycle from any point other than the tip is 
also a state cactus. 

- Consider T>(A, B) : Input Cactus - Defined recursively as 
follows: An input stem with at least one state vertex is an 
input cactus. An input cactus connected to a cycle from any 
point other than the tip is also an input cactus. 

- Consider T>(A) : Chain - A group of disjoint cycles 
(composed by state vertices) connected to each other in a 
sequence, or a single cycle. 

The root and the tip of a stem are also the root and tip of 
the associated cactus. 

Furthermore, recall the following result: 

Theorem 1 ([15]): For an LTI system x = Ax + Bu, the 
following statements are equivalent: 

i) The corresponding structured linear system (A, B) is 
structurally controllable. 

ii) The digraph T>(A, B) is spanned by a disjoint union of 
input cacti. □ 

Note that , by definition, an input cactus may have an input 
vertex linked to several state vertices, i.e., the input vertex 
may connect to the root of a state stem (i.e., input stem) and 
could be linked to one or more states in a chain. 

B. Relation between Maximum Matching and Concepts in 
Structural Systems 

The following results provide a bridge between structural 
systems concepts and graph constructs such as maximum 
matching. These results will be used to characterize the 
minimal dedicated input configurations in Section III. The 
proofs are based on standard graph theoretic properties and 
relegated to the appendix. 

Lemma 1: Let T>(A) = (X,£x,x) and Af* a maximum 
matching associated with B(X,X,£ X<X ), Then, if M* is 
a perfect match, the edges in Af* correspond to a disjoint 
union of cycles in T>(A). □ 

Lemma 2 (Maximum Matching Decomposition): Let 
T>(A) = (X,£x.x) an d M* a maximum matching 
associated with B(X,X,£ x .x)- Suppose Af* consists of a 
non-empty set of right-unmatched vertices. Then, M* (i.e., 
the edges in Af*), together with the set of isolated vertices, 
constitutes a disjoint union of cycles and state stems (with 
roots in the right-unmatched vertices) that span T>(A). 
Moreover, such a decomposition is minimal, i.e., no other 
spanning subgraph decomposition of T>(A) into state stems 
and cycles contains strictly fewer number of state stems. □ 



In other words, the maximum matching problem leads 
to two different kinds of matched edge sequences in Af*; 
sequences of edges in Af * starting in right-unmatched state 
vertices, and the remaining sequences of edges that start and 
end in a matched vertices. These sequences represent state 
stems and cycles, respectively. 

In case a graph is composed of multiple SCCs, we define 

Definition 2: LetZ>(A) = (X,£x,x) and Af* a maximum 
matching associated with B(X, X, £ x>x ). We say that a sub- 
graph of T>(A), denoted by S — (X s , £ x s x s), originates 
a perfect match with respect to Af * if M' — {e : e € 
Af* and e € £ x s iX s} is such that \M'\ = \X S \, i.e, in 
other words, there is a subset of edges Af C Af*, such that 
M' is a maximum matching to B(X S , X s ,£ x s xs). □ 

Definition 3: LetV(A) = (X,£x,x) and Af * a maximum 
matching associated with B(X, X, £x x)- A non-top linked 
SCC is said to be a top assignable SCC if it contains at least 
one right-unmatched vertex (with respect to Af*). □ 

Note that the total number of top assignable SCCs may 
depend on the particular maximum matching Af * (not unique 
in general) under consideration; as such we may define: 

Definition 4: Letf?(^4) = (X, £x,x) and Af* a maximum 
matching associated with B(X, X, £x , x ). The maximum 
assignability index of B(X, X,£x,x) is the maximum num- 
ber of top assignable SCCs that a maximum matching Af* 
may lead to. □ 
III. MAIN RESULTS 

In this section we present the main results of this paper (an 
outline of the proofs is relegated to the appendix), broadly 
centered on the following two issues: 

• Determine the minimum number of dedicated inputs to 
be allocated to ensure structural controllability. 

• Describe the set of all possible minimal feasible input 
configurations (i.e., allocation configurations with the 
minimum number of dedicated inputs) which lead to 
structural controllability. 

The first result relates the minimum number of dedicated 
inputs necessary to ensure structural controllability to the 
structure of the cacti associated with the system digraph. 
Specifically, we have the following: 

Theorem 2: Given the system (state) digraph T>(A), the 
minimum number of dedicated inputs required to ensure 
structural controllability is equal to the minimum number 
of disjoint state cacti that span T>(A). □ 

Theorem [2] reduces the problem of finding the minimum 
number of dedicated inputs to that of finding the minimum 
number of disjoint state cacti spanning T)(A). The next set 
of results are concerned with explicitly characterizing the 
latter number by invoking the relationship (see, for instance, 
Lemma 2) between cacti decompositions and more readily 
computable graph constructs such as maximum matching. 

Minimum Number of Dedicated Inputs 

The following characterization of the minimum number of 
disjoint state cacti spanning T>(A) (and, hence, the minimum 
number of dedicated inputs holds: 



Theorem 3 (Minimum Number of Dedicated Inputs): Let 
T>(A) = (X, £x,x)> where \X\ = n and M* be a maximum 
matching associated with B(X, X, £x ,x)- Moreover, let 
S l = (X l ,£ X i X i), for i = 1, ••• , k, be the non-top linked 
SCCs of the DAG representation of V(A), with X i C X and 
£x i ,x i C £x ,x ■ Then, the minimum number of dedicated 
inputs p is given by 

p = m + j3 — a, (3) 

where 

- m = |V| and V denotes the set of right-unmatched 
vertices with respect to M*\ 

- (3 is the number of non-top linked SCCs; 

-a is the maximum assignability index of 
B{X,X,£ x ,x). " ' □ 

First note that, in Theorem [3] the number of right- 
unmatched vertices (and hence, the minimum number of ded- 
icated inputs p) does not depend on the specific instantiation 
of the maximum matching M* being considered (which is 
not unique in general). Moreover, it may be readily verified 
from the definitions, that if T>(A) is strongly connected, 
we have = l,in Theorem 3, and a may only assume 
two values, or 1, depending on whether m = or 
rn = 1 respectively. As such, Theorem 3 may be simplified 
significantly if V {A) is known to be strongly connected: 

Corollary 1: Let T>(A) = (X,£x,x) be strongly con- 
nected with \X\ = n. Let M* be a maximum matching 
associated with B(X, X, £x,x) and designate by V the set 
of right-unmatched vertices. Then, the number of dedicated 
inputs p 

- is 1 if |V| = (i.e., there is a perfect match), 

- is m if |V| = m > 0. □ 
Characterizations of the required number of minimal ded- 
icated inputs similar to Corollary 1 were presented in [13]; 
however, the strong connectivity of the system digraph is 
necessary for the validity of Corollary 3. In general, when 
T>(A) consists of multiple SCCs, the reduction of Theorem 
3 to Corollary 1 is not valid (as illustrated by the example 
design scenario in Section IV). 

Theorem [3] provides the minimum number of required 
dedicated inputs, hence the minimum number of columns 
in B (each with only one non-zero entry) required to ensure 
structural controllability. We now explicitly characterize all 
such £Ts (up to a permutation of the columns). Each such 
combination will be referred to as a minimal feasible input 
configuration. 

Minimal Feasible Input Configurations 

A minimal feasible input configuration will be denoted 
by a p-tuple < x^ , • • • , Xi > of states, where p corre- 
sponds to the minimal number of dedicated inputs ensuring 
structural controllability (see ([3]) in Theorem 3) and i). £ 
{1, • • • , n} for all k — 1, • • • ,p. In other words, a p-tuple 
< Xi t , ■ ■ ■ , Xi p > corresponds to a minimal feasible input 
configuration if allocating dedicated actuators (inputs) to 
each of the p states Xi k in the tuple leads to structural 
controllability. Also, note that each such p-tuple corresponds 
to a unique canonical B; hence, identifying the set of all 



possible canonical minimal B's is equivalent to identifying 
all such p-tuples of minimal feasible input configurations. 

Also, denote by the set of all possible minimal feasible 
input configurations, i.e., 

={< x i± , - ■■ ,x ip >: x n j£ ■■ ■ ^ x lp and if a 
dedicated input is assigned to each Xi k , where 
ik € {1, • • • , n} and k — 1, • ■ ■ ,p, the resulting 
LTI system is structurally controllable}. 

Note that, by the above definition, a minimal feasible 
input configuration is invariant to any permutation of the 
states in its associated p-tuple representation; a permutation 
of the states in the p-tuple leads to the same dedicated 
input assignment. The following set of results concerns 
the efficient description and enumeration of the set of 
all possible minimal feasible input configurations. The key 
driving factor behind an efficient representation of is the 
existence of subsets J C X, j = 1, ■ ■ ■ ,p, such that is 
almost (to be made precise soon) the Cartesian product of 
the J 's, i.e., ~ 1 x • • • P , i.e., it will be shown that, 
up to permutation and some natural constraints on the J 's, 
a p-tuple < Xi t , ■ ■ ■ , Xi p > belongs to iff . £ Q 3 for 
all j = 1, • • • ,p. Specifically, we have the following: 

Theorem 4 (Naturally Constrained Partitions): Let 
V(A) = (X,£ x ,x), with \X\ = n and M* denoting a 
maximum matching associated with B(X, X, £x ,x)< where 
V = {i>i, • • • , v m } be the set of m right-unmatched vertices. 
Moreover, let S l = (X l ,S x i iX i), for I = {l,--< ,k}, be 
the non-top linked SCCs of the DAG representation of 
V{A), with X 1 C X and £ x i,xi C £x,X- The set of 
minimal feasible input configurations may be characterized 
as follows: 

There exists J C X, j = 1, ••• ,p (where p denotes 
the minimum number of dedicated inputs as in Theorem |3j, 
given by 

- j = 1, • • ■ ,m 

Q 3 ={x : (V - {vj}) U {x} for x £ X is the set of 

right-unmatched vertices for some 
maximum matching of B(X,X,£x,x)}, 

- J "' • •/': (")' U S l , 

ie{i,- ,k} 

such that 

<ii l ,-,ii p >e0 (4) 

if and only if the following hold: 

i) Xi j £ J ij £ {1, • • • , n} for all j = 1, • • • , p and 
i 3 ^ ij, for j ^ ./'. 

ii) Xi j £ J and Xi ., £ 3 for j ^ j' implies that Xi j and 
Xi , are the root to two different minimal state stems, 

3 

iii) for each non-top linked SCC S l , there is at least one 
state variable in < x^,- ■ ■ ,Xi p > that belongs to S'^ 

Note that the sets J are defined on the basis of the specific 
maximum matching M* in consideration. However, as will 
be evident from the proofs, up to a permutation of the 
indices j, j = 1, • • • ,p, the ©■''s are independent of the 
actual instantiation of M* (which may not be unique). We 



refer to the sets J ', for j = 1, • • • ,p endowed with the 
natural constraints as the natural constrained partitions of 
0. Given the description in Theorem |4] we are interested in 
understanding the computational (algorithmic) complexity of 
implementing the natural constrained partitions, as well as 
understanding how to use such characterization to compute 
iteratively a minimal feasible input configuration. This is the 
scope of the next result, stated as follows: 

Theorem 5 (Complexity): Let the hypotheses of Theorem 
4 hold. Then, there exist algorithms of polynomial com- 
plexity (in the number of state vertices) to implement the 
following procedures: 

1) obtaining the minimum number of dedicated inputs; 

2) constructing the natural constrained partitions, 3 's, of 
6; 

3) generating a minimal feasible input configuration itera- 
tively. □ 

First, note that, although by Theorem 5, there exist poly- 
nomial algorithms to constructing the J 's and obtaining 
a minimal feasible input configuration, listing specifically 
all possible minimal configurations may be combinatorial 
(the number of such configurations could be exponential in 
the number of state vertices). The algorithmic procedures to 
obtaining the minimum number of dedicated inputs and the 
natural constrained partitions are provided in Appendix B. 
An example presented in the following section. 

IV. AN ILLUSTRATIVE EXAMPLE 
The following example illustrates the procedure to obtain- 
ing the minimum number of dedicated inputs. Consider a 
set of 6 agents that do share information. Let us assume that 
each agent i has a predefined path y t : K — » M. N , to follow 
parametrized by a time-dependent parameter ji : K — > R, 
such that yi("fi(t)) provides the position of an agent i at 
time instant t. Suppose we are interested in addressing the 
controlled synchronization problem, which consists of some 
predefined parameter specification 7* eR 6 . 

Moreover, suppose that the autonomous system may be 
described as follows: 

71 = -kiji, 72 = -&272, 75 = -^974 , 76 = -^1074, 
73 = -&371 - fc 472 - ^574 , 74 = -^673 ~ &775 ~ ^876, 

(5) 

where fc; 6 E + , i = 1, • • • , 10 are prescribed gains. 

The question that we address is the following: what is the 
minimum subset of agents that we need to assign a dedicated 
control such that any specified 7* is achievable in finite time? 

From (|5j we obtain the structured (autonomous) system 
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The question that we address is the following: what is 
the minimum subset of agents to which we need to assign 
dedicated controls such that any specified is attainable in 
finite time? 

To this end, we consider the following steps: 
Step 1 Compute a maximum matching M* associated with 
B(X, X, Ex ,x)- Figure [ljb) represents in red, the edges 
belonging to the maximum matching M* (note that M* is 
not unique in general). 





where 7 = [71 72 73 74 75 ] T and x denote the non- 
zero entries. System |6} provides the digraph representation 
V(A) = (X, £x,x) as indicated in Fig.[l}a). 
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Fig. 1. a) The digraph T>(A) modeling the inter-dependecy of the 
parameters of the 6 agents; b) The SCCs are depicted by the rectangles 
and the red edges represent a maximum matching (recall Lemma |5J; c) 
Removing the edge in blue from T>(A) forces 71 to be a right-unmatched 
vertex. 



Step 2 Find the minimum number of dedicated inputs 
required to ensure structural controllability. To this end, given 
M* , the set of right-unmatched vertices are V — {73,76}, 
and, hence, in Theorem [5] we have m = 2 and f3 = 2 
corresponding to the non-top linked SCCs. To find a in 
Theorem [5] consider the alternatives for 73, that are in the 
non-top linked SCCs. For instance, let us verify if 71 is a 
possible alternative to 73. To this end, we force -fx to be 
a right-unmatched vertex, by removing all incoming edges 
to 71 on the original digraph (represented in blue in Fig. [2]) 
and compute a new maximum matching. This new maximum 
matching is depicted in Figure[T]-c). Since, the new maximum 
matching consists of the same number of edges, 71 is a 
possible alternative and the SCC containing 71 is assignable. 
Now, consider 76 and let's verify if there exists an alternative 
in the corresponding non-top linked SCC. Recall that 71 
is fixed because it's in an assignable non-top linked SCC 
and our goal is to find the maximum assignability index. 
Let us consider that 72 is also fixed, (by removing its self 
loop), and compute a new maximum matching. This provides 
{7!, 72, 73} as right-unmatched vertices, which implies that 
the maximum matching has one less edge with respect to the 
previous. Since there is only one non-top linked assignable 
SCC, in Theorem [5] we have a — 1. Hence, we need 
p = m + /3 — a = 2 + 2 — 1 = 3 dedicated inputs. 
Step 3 In this step, we characterize all possible feasible 
minimal input configurations (up to permutation) by re- 
sorting to Theorem [4] The generic algorithm (Algorithm 
3) is presented in Appendix B. We start by using Algo- 
rithm 3-(l) (see Appendix B) as a subroutine to compute 
the set of possible alternatives to each right-unmatched 
vertex. Initially set 1 = {71}, O 2 = {76}- In order 
to extend these sets, let's verify which alternatives to 71 
provide maximum matchings with the same number of edges 
as M*. Consider the original B(X, X, £x,x) and fix 7 6 
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Fig. 2. Digraph T>(A) where the SCCs are depicted by rectangles. The 
red edges and state vertices identify the state stems and the blue vertices 
correspond to the dedicated inputs connected to the roots of the state stems. 
In a) we have the feasible minimum input configuration < 71 , 72 , 76 > 
and in b) the feasible minimum input configuration < 72; 71,75 >■ 

(we are just exploring alternatives to 71), by removing 
the edges in Sx.x that end in 76. Moreover, to explore 
if the remaining vertices {72,73,74,75} are viable alter- 
natives, consider for i = 2,3,4,5, B rul6 (X, X, £x.x ~ 
{(.,76) — {(-,7i)}})i ar, d M 7i ' 7e the corresponding (arbi- 
trary) maximum matching. It turns out that |M 72 ' 76 | = 
| M 73,76| _ | M 75,7e| = 4 and |Af73,7e| _ 3^ hence 

O 1 = {71,72,73,75}- Similarly, by fixing 7! we have 9 2 = 
{76,75}. Because /3 — a — 1 in Theorem [3] by Theorem [4] 
we need an extra partition given by 9 3 = {71, 72}- Together 
with the natural constraints we have the characterization of 
all minimal feasible input configurations. 
Step 4 In this step we create iteratively the feasible con- 
figurations from the natural constrained partitions in Step 
3). The natural constraints impose that first we assign a 
dedicated input to at least one state variable to a non-top 
linked assignable SCC. Thus, 

1) Picking 71 from 1 , followed by 72 from 3 leaves 
us with the choice of 75,76 from 2 , leading to the 
minimal feasible configurations < 71 , 72 , 75 > and 
< 7i,72,76 >> respectively. Equivalently, the above 
correspond to following structures of the matrix B (up 
to column permutations): 
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2) Similarly, picking 71 from 3 , followed by 72 from 
1 leaves us with the choice of 75, 76 from 2 , hence 
providing the same configurations as in (1). 

V. CONCLUSIONS AND FURTHER RESEARCH 

In this paper we provided a systematic method with 
polynomial complexity (in the number of the state variables) 
to obtain the minimum number of dedicated inputs (through 
the structural design of B up to column permutation), and 
characterize all possible solutions that ensures structural 
controllability of a given LTI system. By duality, the results 
extend to the corresponding structural observability output 



design. A natural extension of the current framework consists 
of obtaining minimal allocations for general cost constrained 
placement problems, in which actuator-sensor placements 
may incur different costs at different state vertices. This 
problem is more challenging; a natural way to proceed 
is to modify the constructions of the natural constrained 
partitions suitably so as to incorporate the non-homogeneous 
assignment costs. 

APPENDIX A 

Proof of Lemma [7] 

Both sets in the bipartite graph have the same number n of 
vertices, and assume that X = {x , • • • , x n }. Thus, \M* \ = n and 
we may only have the following two possibilities: 

• (a; 1 , a; 1 ) £ M* which is a self-loop (i.e, a cycle) 

• for (x 3 , x k ) G M* there is a sequence of edges such that the 
last edge in M* is given by (.,i J ), otherwise, x J is a right- 
unmatched vertex and we have a contradiction with the fact 
that M* is a perfect match. 

The assertion follows immediately by noting that each of the two 
possibilities above lead to a cycle in T>(A) (by properly sequencing 
the edges in M*) and that the resulting cycles are disjoint in T>(A) 
as the associated edges correspond to a matching. ■ 

Proof of Lemma [2] 

That M* together with the set of isolated digraph vertices span 
T>(A), may be readily verified from the definition of maximum 
matching, see, for example, [15]. 

In the following we prove the minimality of the decomposition 
(through state-stems and cycles) achieved by M* . To this end, note 
that the following generic digraph properties may be verified from 
the definitions: 

(1) Let C denote an arbitrary spanning decomposition of T>(A) 
into disjoint subgraphs of state stems and cycles. Then, the 
edges in C define a matching M c (not necessarily maximum) for 
B{X,X,Sx, x ). 

(2) The root of a state stem belonging to a matching M of 
B(X, X, £x,x) is necessarily a right-unmatched vertex with respect 
to M. 

With the above properties, we now establish the desired minimal- 
ity of M* by contradiction. Assume, on the contrary, there exists 
a spanning decomposition C of T>(A) into disjoint subgraphs of 
state stems and cycles, such that, C consists of strictly fewer state 
stems than M* . Then, by property (1) above, the corresponding 
matching M c consists of fewer state stems than M*. Since, by 
property (2) above, each state-stem in M c corresponds to a unique 
right-unmatched vertex (and similarly for M*), we conclude that 
M c consists of strictly fewer right-unmatched vertices than M*. 
This clearly contradicts the fact, that, M* is a maximum matching. 
Hence, the desired assertion follows. 
Proof of Theorem [2] 

The proof follows by contradiction. Let iV and M denote the 
minimum number of dedicated inputs and state cacti that span 
T>(A), respectively and assume that M 7^ N. Thus, we have two 
possibilities M > N or M < N. 

If M > iV then there is at least one state cacti that, in particular, 
has no dedicated input in its origin. Hence, denoting by B the 
canonical input matrix resulting from the assumed minimal place- 
ment configuration, we note that the augmented digraph T>(A, B) 
is not spanned by a disjoint union of input cacti. Indeed, this 
follows from the fact that a state cactus with a dedicated input 
is in particular an input cactus. 

Now consider the possibility M < N. Similarly, denoting by B 
the (canonical) input matrix, the hypothesis implies that D(A,B) 
is spanned by M disjoint input cacti, which, together with the 
restriction that we have dedicated inputs and the characterization 
in Theorem [T] further imply that TV is not the minimum number of 
dedicated inputs. Thus, the claim follows by contradiction. ■ 



Proof of Theorem [i] 

The proof is lengthy. We sketch the main outline below, the 
intermediate technical details being omitted. 

In order to establish the desired assertion, first note that, by 
Theorem 2 it suffices to show that the minimal number of disjoint 
state cacti that span T>(A) is m + § — a. To this end, let T denote 
an arbitrary decomposition of T>(A) into disjoint state cacti and T 
denote the set of all such T. Clearly, we are interested in obtaining 
min-rgT \T\, where \T\ denotes the number of cacti in T. 

By Lemma 2, each maximum matching of B(X, X, Ex,x) 
generates a T G T and let Tm C T be the subset consisting 
of such T's. Noting the minimality of maximum matching based 
decompositions (again by Lemma 2), it may be argued (we omit 
some techincal details here) that the desired minimation may be 
reduced to that over the set Tm, i.e., p = minreT M T|. 

Noting that each maximum matching consists of m state stems, 
each element T G Tm contains at least m disjoint state cacti. 
Hence, \T\ > m. Hence p > m, i.e., one needs to assign at least m 
dedicated inputs. Now fixing a T G Tm again (and let M* denote 
the associated maximum matching), suppose there are c(M*) non- 
top linked SCCs that originate perfect matches with respect to M* 
and which do not contain right-unmatched vertices of M*. Note 
that, by definition of a, it follows that c(M*) > /3 — a. Clearly, 
for each non-top linked SCC that originate a perfect matches with 
respect to M*, we clearly need to assign an additional dedicated 
input to one of its vertices, because, by definition of non-top linked 
SCC, there is no incoming edge to the states that it consists of. 
Hence, p — mmreT M \T\ = m + rninA/* c(M*) > m + ft — a. 
Now note that, by definition of maximum assignability index, there 
exists a maximum matching (A/*)' such that c((M*)') = ft — a. 
Hence, it follows that p — m + f] ~ a. ■ 

Proof of Theorem |4] 

[=>] First, note that i) holds by definition of €>\ From the proof 
of Theorem [3] it may readily be seen that the natural constraint 
iii) must hold, otherwise the tuple is not a minimal feasible input 
configuration. Secondly, natural constraint ii) holds, otherwise, it 
increases the number of dedicated inputs required to have a minimal 
feasible input configuration. 

[<=] Let < Xi 1 , Xi p > be a p-tuple for which the natural 
constraints hold. Note that, by the hypotheses of Theorem 4, the 
set of p-tuples and the natural constraints are defined w.r.t. a given 
maximum matching M*. From the proof of Theorem [5] we note 
that m entries in the p-tuple correspond to right-unmatched vertices. 
In particular, the natural constraints imply that out of these m 
variables, a are in non-top linked SCCs, and the remaining /3 — a 
arbitrary variables from each of those non-top linked SCCs that do 
not contain any of the previous a right-unmatched vertices. Hence, 
the p-tuple is a minimal feasible input configuration. ■ 

APPENDIX B 

In this appendix we introduce the algorithmic procedures used 
to compute: 1) the minimum number of dedicated inputs; 2) 
the natural constrained partitions, and 3) an iterative method to 
compute a minimal feasible configuration. Using these results 
we then prove Theorem [5] (at the end of this section). In the 
following we introduce two strategies that work as subroutines for 
the main algorithms: STRATEGY A as a subroutine to compute 
the naturally constrained partitions of 6 and the natural constraints 
i)-ii) introduced in Theorem [4j STRATEGY B as a subroutine to 
compute the maximum assignablity index and the natural constraint 
iii). Let B(Si, S2, Es 1 ,s 2 ) be a bipartite graph and M* an asso- 
ciated maximum matching M* (computed, for instance, using the 
Hopcroft-Karp algorithm [14]). Now, we show how to compute a 
maximum matching such that a specific vertex is right-unmatched 
(STRATEGY A), or, how to compute a maximum matching with a 
specific matched edge in it (STRATEGY B). 



ALGORITHM 1: Find the right-unmatched vertices that 
belong to non-top linked SCCs 

Input: Set of right-unmatched vertices V, B(X, X,£x x) and 

S\ i = l,.../9. 
Output: Set V of right-unmatched vertices in the non-top 
linked SCCs 

V = {}; 

for each z G V do 

for each v G S 1 U ■ ■ ■ U S ~ V do 

Compute a maximum matching My associated with 
B(X,X,£x,x — {(x, v) : x G X} — {(x, w) : x G 
X,w G V- {z}}); 
if |M„| = m then 
I V = V U {v} and EXIT FOR 



STRATEGY A - Find a maximum matching with some vertex 
as a right-unmatched vertex 

By definition, a right-unmatched vertex is a vertex that does 
not belong to the tip of a directed edge from S\ to S2. To this 
end, let B(Si, S2, Es 1 ,s 2 ~ {( s i, v ) '■ s i G Si})be the graph 
(bipartite) obtained by removing all the edges ending in v; then 
any maximum matching M v of the latter graph contains u as a 
right-unmatched vertex. Note that in our case (recall Lemma [2jl, 
each minimal state stem has its root in a right-unmatched vertex, 
hence we are interested in the cases where \M V \ = \M*\, in other 
words, when the number of state stems does not increase. 

STRATEGY B - Find a maximum matching with a given 
specific edge being a matched edge 

To find a maximum matching for which a prescribed edge 
e = (ei,e2) is a matched edge, let M e = {e} U M', in 
which M' corresponds to a maximum matching associated with 
B(5i-{ei}, 5 , 2 -{e 2 }, E Sl ,s 2 -{(ei, «a) : "2 G S 2 }-{(vi, e 2 ) : 
vi G Si}). It may be seen that M c is in fact a maximum matching 
and contains e. o 

Finding the minimum number of dedicated inputs 

Recall Theorem [5] First compute a maximum matching M* 
associated with B(X, X, £x,x)', we then obtain a set of m right- 
unmatched vertices V = {«!,••• ,v m }. Secondly, compute the 
DAG of T>(A) and denote by «S% i = 1, • • ■ ,P the non-top linked 
SCCs (and f3 the same as in Theorem [3](. Now, sequentially try to 
find a maximum matching such that the right-unmatched vertices 
are in the non-top linked SCCs and keep record to which non- 
top linked SCC they can be assigned. Based on this record, the 
maximum assignability index (a in Theorem [3} may be obtained 
as shown in the following. 

Algorithm 1 provides us with the right-unmatched vertices that 
can be in the non-top linked SCCs. Now we need to maxi- 
mize the number of right-unmatched vertices in different non-top 
linked SCCs. This procedure, in turn, will be used to provide 
us with the maximum assignability index. To this end, consider 
B(I v ,Is, Ei v , i s ) where 

• Iv = {1, •• ■ >|^|} is the set of indices of the vertices in 
V — {«!,••• ,t!|yi}, i.e., the right-unmatched vertices that 
are in the non-top linked SCCs (provided by Algorithm 1); 

• Is = {1, • • • , 0} denotes the indices of the non-top linked 
SCCs; 

• Ej v j s is the set of edges e = i G Iv, j G Is that 
codifies the possibility of a vertex Vi belonging to a non-top 
linked SCCs S J (provided by Algorithm 2). 

By definition, the maximum assignability index, may now be 
readily computed by obtaining a maximum matching associated 
with B(Iv, Is, Ei Vi i s ). Finally, the minimum number of dedicated 
inputs is obtained asp = m + /3 — a by Theorem 3. 



ALGORITHM 2: Determine E IvJs 

Input: V from Algorithm 1 and B(X, X,£x x) and 

Sf j = !,■■■ ,P 
Output: Ei v ,i s (possible assignments between dedicated 
input and non-top linked SCCs) 

Ei v j s = {}; 

for each v € V = {vi, ■ ■ ■ , U|vi} do 

for each w e S 1 U ■ • ■ U S fi - (V - {v}) do 

Compute a maximum matching M associated with 

B(X,X,£ x .x - {(x,w) : x £ X} - {(x,vi) : x £ 

X, Vl e (V — {«})}); 
if |M| = m then 
(consider v = Vi) 
for j = l,- ■■ ,<3 do 
if iu £ 5 J then 
I Sj Vl i s =S Jv , Js U{(t,j)} 



Generating a minimal feasible input configuration 

Given a maximum matching M* of T>(A), in this section we 
provide an iterative algorithm to generate a minimal feasible input 
configuration, which is a tuple < Xi x , • • • , Xi p > invariant to 
permutation. Broadly, the idea is to iteratively choose a variable 
from the natural constrained partitions (note that we may select 
the first variable from any of the partitions) and update the natural 
constraints based on which a new variable may be chosen. 

Specifically, the strategy to create a minimal feasible input 
configuration consists of the following: first, chose the variables that 
are in the non-top linked SCCs that ensure the maximum number 
of top assignable SCCs. Secondly, select one state variable for each 
non-top linked SCC that was not selected in step one. Finally, the 
third step consists of selecting the variables corresponding to right- 
unmatched vertices with respect to M* , that were not used in step 
one. 

Proof of Theorem [5] 

1) Minimum number of dedicated inputs: It consists of computing 
an initial maximum matching in 0(^/|A'j|£;f.;r|) [14] to find 
m. To determine the DAG (may be implemented by executing 
Depth-First Search twice, with complexity 0(1^" + jf^A-j) [14]). 
The maximum assignability index is computed by executing Al- 
gorithm 1 and Algorithm 2; the first algorithm has complexity 
■mNO(^/\X\\£x,x\), with N (< N) being the number of state 
vertices in the non-top linked SCCs. Algorithm 2 has complexity 
0(iV| V|) with \V\ being at most m. 

3) Step (1) (Algorithm 3) executes the maximum matching at most 
ap\X\ times, hence its complexity is C(ap|< ; f|-\/|A'||£'A',A'|). Step 
(2) (Algorithm 3) uses at most p\X\N executions of the Algorithm 
2. Step (3) (Algorithm 3) consists of a choice of a single state 
variable by the designer, hence of constant complexity. Step (4) 
(Algorithm 3) consists of j3 — a executions of a constant complexity 
procedure, hence, is of constant complexity again. Finally, Step (5) 
(Algorithm 3) executes p — j3 times a procedure of complexity 
0(p\X\<J\X\\£ x ,x\). 

2) The natural constrained partitions can be computed by using 
single execution of Algorithm 3. ■ 
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